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The goal of this paper is twofold: to explore the response of classical charges to electromagnetic 
force at the level of unity in natural units and to establish a criterion that determines physical 
parameters for which the related radiation-reaction effects are detectable. In pursuit of this goal, 
the Landau-Lifshitz equation is solved analytically for an arbitrary (transverse) electromagnetic 
pulse. A comparative study of the radiation emission of an electron in a linearly polarized pulse 
for the Landau-Lifshitz equation and for the Lorentz force equation reveals the radiation-reaction- 
dominated regime, in which radiation-reaction effects overcome the influence of the external fields. 
The case of a relativistic electron that is slowed down by a counterpropagating electromagnetic wave 
is studied in detail. We further show that when the electron experiences acceleration of order unity, 
the dynamics of the Lorentz force equation, the Landau-Lifshitz equation and the Lorentz- Abraham- 
Dirac equation all result in different radiation emission that could be distinguished in experiment. 
Finally, our analytic and numerical results are compared with those appearing in the literature. 
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PACS numbers: 03.50.De,33.20.Xx,41.60.-m,41.75.Jv 

I. INTRODUCTION 

A. The radiation-reaction force and 
Landau-Lifshitz equation 

The equation of motion for an electron with charge — e 
and mass m in an external electromagnetic field is given 
by the Lorentz force (LiF) equation 



mil" = -eF°"\ip, 

where F"'' is the electromagnetic tensor, u" 
is the four- velocity of the charge, 7 = 1/vT" 



(1) 

jy_(l,t/) 
u^ and 
the dot represents differentiation with respect to proper 

time T. We use the metric convention (H ) and 

units in which the speed of light is unity c = 1. Maxwell 
equations together with the LF equation, imply that the 
rate at which energy is emitted by an accelerated charge 
with respect to laboratory time is [l[ 



TZ 



dE 
'dt 



-e^u^iin 



(2) 



lab 



The right-hand-side of Eq. ([2]) is a positive-definite 
Lorentz invariant that vanishes if and only if fi'^ = 0, 
and therefore an accelerated charge emits radiation and 
loses energy relative to any Lorentz observer. 

The inertial reaction due to the energy-momentum loss 
exhibited by Eq. ([2]) is not accounted for in the LF equa- 
tion 1^. Consequently, the Lorentz force equation with 
a prescribed electromagnetic field is only an approxi- 
mated description of the electron motion limited to cases 
in which the radiation emission is small. The precise 
meaning of "small" will be established in Sec. IIVI where 
the radiation-reaction (HRJ dominated regime criterion 



suggested by analysis [1, 3 is verified by numerical eval- 
uation of the radiation. 

During the development of the Lorentz-Maxwell the- 
ory of electromagnetism, there has been a long search for 
an improved classical equation that comprehensively de- 
scribes the motion of a radiating charge in a prescribed 
electromagnetic field. The most prominent equation that 
was suggested after the introduction of quantum theory 
is the Lorentz-Abraham-Dirac (LAD) equation Q 



mu" = -eF'^Pup + ruTo [11°" + u^m"] , 



where 



"^o "" q — 3 



(3) 



(4) 



is a constant with dimensions of time, whose numeri- 
cal value for the electron is tq = 6.24 x 10~^^s. Other 
models have been introduced by Eliezer [5| , Landau and 
Lifshitz [3|, Mo and Papas [0], Caldirola Q, Yaghjian Q 
and Sokolov et al. ^^. 

Of particular interest is the Landau-Lifshitz (Lh) equa- 
tion (also known as the reduced LAD equation) first pre- 
sented in [3[ 



-ei^"%-ero n^U;3M 



.1. [F'^Ppfs^u^ - F^^F^su^Ufiu''] |. 



(5) 



The LL equation is often selected for further study be- 
cause it is the only equation from the above list that is 
equivalent to the LAD equation up to first order in tq [12| , 
avoids the nonphysical solutions of the LAD [l3| and has 
known analytic solutions |14l4l6l |. A rigorous derivation 



of Eq. ([5]) using perturbation theory was recently given 
from considerations of energy and momentum conserva- 
tion [l3|. 



B. Line of approach and objectives 

In this paper we explore the dynamics of a charged par- 
ticle exposed to an ultraintense pulsed laser field. The 
objective will be to understand how RR impacts the nor- 
mal Lorentz dynamics and to derive criteria which delin- 
eate the domain of validity of the LF. For a thorough dis- 
cussion of the electron dynamics and radiation emission 
in laser fields with a small radiation reaction correction, 
the reader should consult Ref. [3- The assumption of 
small radiation-reaction breaks down as the acceleration 
reaches unity in natural units. For acceleration unity, the 
familiar Lorentz force dynamics cannot hold anymore, 
and RR effects dominate the motion of the charge as will 
be discussed in depth in the body of this work. 

The investigations here are intended to help identify 
the onset of new dynamics beyond that found in classi- 
cal electromagnetism, and not so much to test which of 
the LF generalizations is more accurate. Moreover, since 
the theory of quantum electrodynamics is founded upon 
classical electromagnetism, the completeness and depth 
of our understanding of charged particles under high ac- 
celerations remains both qualitatively and quantitatively 
uncertain [l9|. In particular, dynamics at the critical 
electromagnetic field strength Ec = m?C' /eh — both clas- 
sical and quantum — should be reconsidered [201 . At this 
limit quantum effects are believed important |2ll . |22J [. 
though wc withhold judgment on the matter and con- 
sider a thorough analysis of classical predictions a rele- 
vant basis for further exploration. 

To reach acceleration unity with current laser systems, 
we boost the intensity of the laser wave by colliding it 
head-on with relativistic electrons (Fig. [IJ. In the in- 
stantaneous rest frame of the electron the laser fields are 
greatly enhanced, as was demonstrated numerically in 
[J, [23]. Since the LL equation is considered to be the best 
available candidate to account for the radiation-reaction 
effects, we solve the LL equation and analyze the dynam- 




FIG. 1: The experimental setup, represented qualitatively. 
The radiation propagates mostly in the direction of the elec- 
tron and is projected onto the screen on the left. 



ics according to this model. 

The LL Eq. ([5]) is nonlinear in the electromagnetic field 
tensor F"^ and in the four- velocity u" . Because of this 
nonlinearity, it has resisted for a long time an anal ytic 
solution, except for simple cases: constant magnetic |24l . 
l25j and constant magnetic plus electric fields for which 
the LL equation reproduces the dynamics of the LAD 
equation [26| , a circular orbit [ij] and the nonrclativistic 
motion in a Coulomb potential ^M,- Recently ^M,^ a 
closed-form solution for LL was obtained for a plane wave 
and studied in applications [27 1. 

In this paper we solve for the motion of a charged par- 
ticle in a transverse wave: a monochromatic electromag- 
netic wave for which the wave fronts (surfaces of con- 
stant phase) are infinite parallel planes of arbitrary am- 
plitude. This includes many useful physical scenarios as 
special cases, such as constant crossed electromagnetic 
fields, linearly polarized and circularly polarized plane 
waves, and an electromagnetic pulse in space and time 
with arbitrary shape. The solution here is obtained by a 
method independent of that seen in |16| , and we provide 
additionally an analysis of the radiation emission and its 
angular distribution. Analytic and numerical compar- 
isons of our solution with those of [la, 123 ^'"^ provided 
in Appendix [Cl 

In Sec. In] we discuss the physical meaning and math- 
ematical properties of the transverse wave. Section IIIII 
includes the derivation of the analytic solution of the LL 
equation for this case. In Sec. IIVI wc compare the radia- 
tion emission with or without the RR force, and discuss 
in what conditions the LL equation can be probed ex- 
perimentally. In Sec.|V]we study our two main examples 
of linearly polarized and circularly polarized plane waves 
and compare them to the known solution of the LF equa- 
tion. 



II. THE TRANSVERSE WAVE 

We consider the motion of a charged particle in a trans- 
verse wave: a monochromatic electromagnetic wave for 
which the wave fronts arc infinite parallel planes of arbi- 
trary amplitude. The four-potential of such an electro- 
magnetic wave is 



A"(x)-AoRe[e"/(C)], 



(6) 



where Aq is the maximal amplitude of the wave, e" is the 

(complex) polarization four- vector, fc" = fi^,^") is the 

propagation four- vector and / is an arbitrary (complex) 
function of ^ = fc^a;^ = fc • a: that represents the shape of 
the wave. 

To keep the dependence on the intensity of the elec- 
tromagnetic wave explicit, we introduce the normalized 
four-potential 



i"(x)=Re[e"/(0] 



(7) 



so that the four-potential is A" = AqA". This will prove 
useful in the analysis of the radiation emission in Sec. lIVI 
We take a polarization four- vector e" and a propaga- 
tion four-vector fc" that satisfy 



= 0, 
= -1, 



(8) 



and the transverse condition 



0. 



(9) 



This means that the polarization of the wave is orthog- 
onal to the direction of wave propagation. Namely, the 
electric and magnetic fields are perpendicular to the di- 
rection of energy transfer. 

The transverse condition Eq. © guarantees that the 
electromagnetic potential satisfies the Lorenz gauge con- 
dition 

daA°' = 0. (10) 

The electromagnetic field tensor is 

F^'^ix) = d"A'^-d^A°' (11) 

= AoRe [{k^e^ - k^e") f (0] 
= k^'A"^ -k'^A'"', 

where the prime denotes differentiation with respect to 
the variable ^. Equation ([5]) and the transverse condition 
Eq. ^ imply that the field tensor satisfies the following 
identities 



tion. In this case 

/(C) = 
A = 



(0,1,0,0), 

(^,0,0,fc), 

sin(^-^o), 

—Aq sin(fcz ^ Lot + ^q)x^ 



(13) 



E = —LuAQCOs{kz — iot + ^Q)x, 
B = —kAo cos{kz — ujt + ^o)y, 

where ^o is the phase of the wave. The choice of phase has 
an important physical significance, as it determines the 
intensity and the direction of the electromagnetic wave 
as it initially hits the particle. 

b) Circularly polarized plane wave propagating in 
the positive z direction with positive helicity. In this case 



1 



(0,1,-1,0), 



(14) 



fiO 

A 
E 



V2 
(w,0,0,fc), 

Aq [cos(fcz — wi + ^0)2; — sin(fcz — ujt + Co)y] , 
—u)Aq [sin(A:z — ujt + ^q)x 
+ cos{kz — ujt + £,o)y] , 
B = —kAo[~cos{kz — ojt + £^o)x 
+ sin(fc2 — ujt + £,o)y] . 

This example of a circularly polarized plane wave also 
demonstrates the importance of allowing e" and /(^) to 
be complex. 

III. THE SOLUTION 



fc„i^"^ 


= 0, 


(12a) 


fcaF.f 


= 0, 


(12b) 


e.F"^ 


= -(e.A')fc^ 


(12c) 


£.^.f 


^ ^{e-A")k^k„ 


(12d) 


F'^^u, 


= (wA')fc"-(fc-w)A'", 


(12e) 


^'^^Fp, 


= -fc"fc^(A')'. 


(12f) 



Equation (|12ap reveals that the potential ([6|) generates a 
very special field configuration, in which the wave four- 
vector k" is orthogonal to each of the row/column four- 
vectors of the Faraday field tensor F"^ . We will see later 
that this implies that the four-acceleration u°' is orthog- 
onal to fc" in the case of the LF Eq. ([1|). However, the 
orthogonality of the four-acceleration and the wave four- 
vector will cease to hold once RR terms are included, in 
the case of the LL Eq. ([5|). 

Two concrete examples of the transverse wave Eq. 1^ 
we will address in Sec. FVlare 

a) Linearly polarized plane wave propagating in 
the positive z direction with polarization in the x direc- 



The gist of the method used to solve the LL equation 
is to introduce a change of variables similar to the one 
used in [g^ for the LF Eq. ([T]) . We use the phase S, = k-x 
instead of the proper time r as the independent variable. 
This means that ^ is related to the proper time t by the 
relation 



-— = k ■ u 
dr 



(15) 



since the wave vector fc" is fixed. Writing the LL Eq. ^ 
in terms of ^, one obtains 



(fc • u) u" 



_jipap 



m 
e 



up - -rJFfupu^ 



[F'^^Fp^u' - upF^-'F^su'u'^]], (16) 



where the prime denotes differentiation with respect to 
the new variable ^. Equation (|16p is a differential equa- 
tion in u°' that contains the terms e ■ u and A: • w as we 
can sec from identities ([T2l) . Contracting Ec^. (fT6|) with 
ka and using the field identities (|12|) gives 



{k-u)' = T^al{k-uY{A')\ 



(17) 



where 



and dividing by (fc • uf' we have 



ao 



cAq 



(18) 



is a positive dimensionless constant measuring the inten- 
sity of the electromagnetic wave. 

Dividing by {k ■ u)^ and integrating, we have 



k- 



k ■ uq 



1 -roag(fc--uo)V'(C)' 
where we defined the first structure integral 



^(0 = / A'{y) 



dy. 



(19) 



(20) 



Relation ((TS|) allows integration of Eq. P^ . obtaining 
an explicit expression for r as a function of ^, namely 



r(0 



e 



k ■ Uf) 



Jo 



(21) 



If A'" is a spacelike vector (e.g., when the time com- 
ponent of the polarization four- vector vanishes, e'^ — 0), 
then when ^ > the function ijj is a non-negative func- 
tion. Thanks to Eq. dUl) we now sec that H > 0. There- 
fore ^ is an (increasing) monotone function of r, and the 
change of variables r — s- ^ can indeed be used for r > 
(notice, however, that it might run into a singularity if 
the proper time is negative). 

There is a caveat in Eq. (fT9)) that the reader should 
be aware of. Since we perform the integration with re- 
spect to ^ and not t, the constant of integration should 
be determined by setting ^ = 0. Given uq = u{^ = 0)i in 
general Ug is not the initial four-velocity of the particle. 
In order to remedy this situation and to minimize confu- 
sion, we now choose the coordinate system such that at 
T = 0, the particle is at the origin and therefore ^ = as 
well. Since the change of variables was proved to be one- 
to-one, this guarantees that uq = u{^ = 0) = u{t = 0), 
so for this particular coordinate system uq is indeed the 
initial four- velocity of the particle. 

We continue by contracting Eq. P^ with Sa and using 
the field identities P^ in a similar fashion, which gives 

(fc • u){e ■ uY = ao{k ■ u){e ■ A') + Toao{k ■ u)^{e ■ i") 

+Toal{e-u){k-u)^{A'f. (22) 

This is a nonhomogeneous linear differential equation for 
e ■ u which can be solved given k ■ u. The last term 
in Eq. ([22|) is (e ■ u){k ■ u)' as we can see from Eq. (fTT]) . 
Therefore, one can write 

{k-u){e-u)' - {k-uy{e-u) ^ (23) 

ao{k ■ u){e ■ A') + Toao{k ■ u)^(e • A"), 



e ■ u 



ao- h Tofloe • A . 



■.k ■ u/ k ■ u 

Substituting Eq. (J19p for k ■ u and integrating yields 



(24) 



e • u 
k • u 



£•"0 , e-{A~ Aq) 3 

-I- ao ; Toage • X (25 

k ■ Uq k ■ Uq 

+Toaoe-{A' -A'q), 



where here again Aq — A°'{0) and we defined as a second 
structure integral the four-vector 



x"{0 



A'"(y)^(y)dy. 



(26) 



Multiplying Eq. ([25|) by /' and /" and taking the real 
part gives us 

A'-u A' -up , A'-{A-Ao) 

-, = -; hao (27 

k ■ U k ■ Uq k ■ Uq 

-To4{A' ■ x) + Toaoi' • {A' - A'„) 



and 



A"-u A" -up A"-iA^Ao) 

—, = —, r ao ; [2^) 

k ■ U k ■ Uq k ■ Uq 

-roag(i"-x) + Toaoi" -(i'-io)- 

We are finally ready to integrate the LL Eq. (|16l) . 
Writing it explicitly in terms of the four-potential using 
Eq. (fTTj) gives 



{k ■ u)u" 



-ao 



-Toflo 



-t-oOq 



{A' ■ u)k°' - (fc • u)A'°' (29) 

{k ■ u){A" ■ u)k°' ^ {k ■ uf A""" 
[k ■ u)A''^k°' - [k ■ ufA''^u°' 



Similarly to what we had earlier, Eq. ((T7)) reveals that 
the last term is [k ■ uyu°', therefore 



k ■ u 



ao 



k ■ u 



-A" 



Toao^" 



ao 



A'-u 
[k ■ uY 



A" -u 



Toflo- 



TqUq 



k- 



fc",(30) 



where we separated terms that are parallel to e" in the 
first brackets from terms that are parallel to fc" in the 
second brackets. Substituting Eqs. dH]), (gT]) and (pgj) . 
integrating and collecting powers of tq, we finally obtain 



k ■ u 
k ■ uq 



aoiA"~A^) 



k ■ Uq 



k" 
k ■ Uq 



aQ{A- Aq) -UQ + aQ 



,{A~AQf 



3,, a 



+ ik-u)TQ\ ao(^'" - ^o") - aoX' 



+k''{k-u)Tn -al 



:iA'-A'Qy 



^2 



,6% 



+ aUA'^A'Q)-x + at^-4^\. 



-ao(i' - A'q) ■ uq - al^/j ~ al{A - Aq) ■ {A' - A'q) + agx • uq + 4(A - Aq) ■ x 



(31) 



r 



This is an analytic expression for the four- velocity u" 
as a function of the variable £,. It is verified by direct 
computation that u^Ua = 1 as expected. Moreover, it is 
given in a manifestly covariant form that is valid for any 
reference frame. The presence of the four-potential A" 
makes the solution not manifestly gauge invariant. Nev- 
ertheless, the solution is invariant under gauge transfor- 
mations satisfying the transverse condition Eq. 0, i.e. 
transformations of the form A°'{^) -^ A°'{^) + k°'A{Q. 

Equation (f3T|) together with Eqs. (fT9|) and (|2T|) provide 
a complete description of the velocity of the particle as a 
function of its proper time. Once the structure integrals 
Eqs. (f20| and ([26| are evaluated, Eq. ([3T|) forms an an- 
alytic solution of the LL Eq. ^ for the transverse wave 
potential Eq. ([5]) in a closed form. The equivalence of the 
solution presented here with [16| is proven analytically 
in Appendix [C] A direct comparison of numerical results 
for the physical situation studied in [231 i^ also presented 
there. 

Since the LF Eq. ^ is obtained from the LL Eq. ([5]) in 
the limit tq — )■ and the field is continuous, the solution 
of the LL equation contains the solution of the LF Eq. ([T]) 
as a special case. In the limit of tq — )• 0, we see that the 
solution Eq. (|3ip reduces to 



+ ao{A''~A^) 



(32) 



k° 



k ■ Uq 



aQ{A - Aq) ■ uq + , 



■.{A-AQf 



which is the known solution [1^ to the LF Eq. ^. In 
this case £, and r have linear dependence, and Eq. ((2T|) 
reduces to 



^(r) = (fc-uo)T. 



(33) 



Comparing the solution Eq. ((3T|) to the LL equation 
with the solution Eq. ([5^ to the LF equation we see that 
To and Tq appear in the former. Terms that are quadratic 
in To are in the direction of the wave vector fc". The linear 
terms in tq were separated intentionally in Eq. (|3ip to 
terms in the direction of e" (last term in the first line) and 
terms in the direction of k" (whole second line). Most 
of the new terms that produce the deviation from the 
LF prediction are in the direction of propagation of the 
wave, as one would expect. 

Despite the appearance of (at most) quadratic terms 
in the solution Eq. pi[) . it is important to realize that the 



dependence of m" on To is more complicated than what 
meets the eye. The projection k-u can be expanded in an 
(infinite) Taylor series in tq using Eq. ([T5)) . from which 
it is seen that u" contains an infinite number of powers 
of To. This all-orders expansion does not imply that the 
LL equation gives an exact description of the radiating 
charged particle. The LL Eq. (j31[) is an approximation 
taking into account the radiative energy loss of the par- 
ticle only from the zeroth order (Lorentz) motion given 
in Eq. ([1]) . For a more complete discussion of the physics 
and assumptions involved in this expansion, see [17| . 

The solution Eq. pip to the LL equation contains pow- 
ers in the intensity up to Uq, four powers higher than Qq, 
the highest power appearing in the solution Eq. (|32p to 
the LF equation. Consequently, increasing the intensity 
of the electromagnetic wave brings the dynamics closer 
to the RR dominated regime, in which the predictions 
of the LL model differ significantly from those of the LF 
equation as we now discuss. 



IV. RADIATION EMISSION 

A. Radiation reaction 

The rate of energy-momentum loss due to radiation is 
normally computed by the Abraham-Heaviside formula 






2)w". 



(34) 



However, a careful examination of its derivation [l| shows 
that Eq. (p4| depends on the equation of motion, which is 
assumed to be the LF Eq. ([T|) . This renders a serious diffi- 
culty in the computation of radiation emission for the LL 
equation, as it currently does not have a corresponding 
expression for the rate of energy lost in radiation. The 
difficulty stems from the implicit description of radiation 
reaction in terms of the Maxwell equations and thus ab- 
sence of a complete theory of chargcd-particle dynamics, 
as would be encoded in a single unified action principle. 
We are therefore compelled to proceed by retaining the 
usual expressions for radiation emission obtained from 
the Maxwell-Lorentz theory. The present work should 
be considered as a case study for the radiation emission 
where RR effects are present, computed within the limi- 
tations of the current theoretical framework, as no other 



approach yet exists. 

The rate of energy-momentum loss Eq. ([M]) is a four- 
vector that requires the knowledge of the four-velocity 
and four-acceleration of the particle. The four-velocity 
was already given in Eq. pip , and the four- acceleration 
is supplied in Appendix [BJ Since the four- velocity and 
four- acceleration arc different in the Lorcntz and in the 
LL dynamics, the radiation emission produces different 
results as well. The rate (with respect to proper time) 
at which energy is emitted from the particle is just the 
0-component of the four- vector in Eq. ([M)) . In order to 
obtain TZ, the rate at which energy is emitted with respect 
to laboratory time, we notice that 



n 



"dT 







drdP 
dt dr 



' 2-2 

-e u , 



(35) 



where wc used Eq. (|34|) in the last equality. This is noth- 
ing other than Eq. ([2]) which was stated earlier in the 
introduction. Evaluating TZ using the four-acceleration 
Eq. (JB1[) yields the rate at which energy is radiated from 
the charge with respect to the laboratory time according 
to the LL equation 



n 



r 



{k ■ uY 



{agi" 



+Toik ■ uo) halA" ■ A' - 2af^4'A''^ 



+T^{k-Uof 



A< A'2\2 



ha^(A 



alA"^ - 2a^'ipA' ■ A 
-2 



-0^ zao'( 

k ■ u f k ■ u 



k ■ uq \k ■ uq 
+T^{k-u){k-u„f{A'yali'(2~2 

+T*{k-ufik-uofiA'^ral^^]. 



k ■ u 



(36) 



With the aid of Eq. ([T^ , it can be written as 



2 2 (fc • ")' 



3 {k-uoY 

oqA' + [k ■ wo)ToaoA" — {k ■ uo)T'oaoV'^' 



(37) 

n 2 



-{k-u,YT^at{A'') 



Since the four-vector in the brackets is a spacelike four- 
vector, this is a manifestly positive-definite expression 
[this is also evident from its definition in Eq. psp ]. By 
taking the limit tq — >■ one obtains the radiation rate in 
the absence of RR effects 



Tl^--e\k-u,falA'\ 



(38) 



This is the rate at which energy is radiated from the 
particle in the case of motion according to the LF Eq. ([Ij . 
Comparing Eq. (PSI) with Eq. ([55]) allows us to determine 



the condition in which the predictions of the LL model 
depart from the predictions of the LF equation. The 
particle has an initial velocity of Wq = 7o(l, uq), and 



k-uo = 7o(w- k ■v'o). 



(39) 



If the wave vector k and the initial 3-vclocity vq point 
in the same direction, k ■ uq < uj. However, when the 
wave vector k and the initial 3-velocity vq are in oppo- 
site directions, k ■ uq has order of magnitude 0170. This 
corresponds to a head-on collision between the electron 
and the laser beam, which thanks to the factor of 70 en- 
hances the radiation-reaction effects greatly as we will see 
next. When oq ^ 1, the leading order correction to the 
LF equation will be the tq term, in which the dominant 
term is that proportional to Qq. Wc therefore enter the 
RR dominated regime when the corrections supplied by 
the LL equation become of the same order of magnitude 
as the ones in Eq. (pS)) . namely as 



al 



(wro)7oao. 



(40) 



For an electron traveling into a laser beam with wave- 
length A ~' 942 nm (or equivalcntly frequency u) ^ 2 fs~ ) 
we have lotq ~ 10~®, and we enter the RR donfinatcd 
regime as 



ToOq 



10« 



(41) 



Obtained on the basis of an analytical solution, these 
conditionsprovide nontrivial validation of the conditions 
found by [3| and [J]. Furthermore, with an explicit solu- 
tion in hand we can go further and verify the condition 
(j40p by directly comparing experimental observablcs pre- 
dicted by the LF and LL equations. 

We sec that there is an inverse squared relation be- 
tween the initial energy of the particle and the intensity 
of the laser. Notice, however, that once we are deep in- 
side the radiation dominated regime, i.e. when 



Qq < (wTo)7oao 



(42) 



higher order terms in tq begin to dominate the emission 
from the particle and not the term linear in tq. This 
still results in a different radiation emission than the one 
predicted by the LF equation. 

Figure [5] demonstrates the validity of Eq. (|¥T|) in the 
case of a head-on collision between an electron and a 
linearly polarized plane wave. The density at each value 
of the intensity oq and 70 was computed numerically over 
one period of the wave (in ^) by considering the statistical 
deviation in the energy of the particle predicted by the 
LL equation and the LF equation. It is a value between 
and 1 given by 
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The result is striking, as Fig. [5] shows that already with 
current laser systems RR effects can be tested. In fact, 
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FIG. 2: A density plot of Eq. (|43p for a LP wave with wave- 
length A = 942 nm. Above the dashed (blue) line given by 
Eq. (|47p . the predictions of the LL equation differ from those 
of the Lorentz-Abraham-Dirac equation. Along the solid 
green line the electron experiences a force of order unity [m^] 
[see Eq. (fig)) ]. The Schwinger critical field Ec = rr?(? jeh is 
indicated in the bottom right corner. The solid circles are the 
values that are studied in Sec. IVl 



the figure shows that RR effects can already be tested 
for high-intensity lasers with oq = 100 and 70 = 1000 (a 
"modest" initial energy of i?o = 0.511 GeV), as is demon- 
strated in the next section. The boundary of the shaded 
domain is in agreement with condition Eq. (j4ip . 

Recall that the LL Eq. ([5]) was derived from the LAD 
Eq. dSl) using perturbation theory Q ■ The earlier edition 
[2g | makes it clear that the radiation-reaction forces of 
the LAD and LL equations agree (perturbatively) when 
the Lorentz force dominates. However, solutions to the 
LAD or LL equation need not be equivalent in general, 
though a relationship can be established in the case that 
the radiation-reaction force remains subdominant |13| . 

Specifically, LAD and LL predictions can disagree 
when radiation-reaction forces dominate the Lorentz 
force. In fact, a difference is possible in a region of param- 
eter space where classical dynamics are valid, as is seen 
by comparing the dominant term in the LL force Eq. ([5]) 
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FIG. 3: The Lorentz invariant acceleration y/—u°''u,a in nat- 
ural units for a circularly polarized laser wave with ao = 100 
and initial 70 — 1000. The solid red line is the acceleration in 
the LF case, while the dashed blue line gives the acceleration 
according to the LL. 



For wavelengths A ~ 942 nm we have ojtq ^ 10^^, and 
the LL equation is compatible with the LAD equation 
when 



ao7o < 10^ 



(47) 



As noted also in Sec. 76 of Q, this condition differs from 
Eq. (|41[) : In the former, oq is proportional to the square of 
the inverse of 70 and in the latter 70 is proportional to the 
square of the inverse of ag. Consequently, there are values 
of ao and 70 for which radiation reaction dominates the 
dynamics. 

We can therefore arrange experiments for which the 
LF Eq. HD, the LAD Eq. (P and the LL Eq. © may all 
predict different particle dynamics. For these values we 
cannot only measure RR effects, but can also distinguish 
the LL equation from the LAD equation. The dashed 
(blue) line in Fig.[2]represents the quantitative condition 
in Eq. (|^7)l . for values at which the RR force is 1% the 
magnitude of the LF, and above it, LL and LAD dy- 
namics will be distinguishable. For example, the point 
marked at intensity oq = 100, initial 70 ~' 10'^ suffices to 
probe the RR dominated regime in which the LL Eq. ^ 
is incompatible with the LAD Eq. ([3]). 

Since the RR force stems from radiation emission, and 
radiation is emitted only if the particle is being accel- 
erated [see Eq. ([2])], the real key physical property that 
is responsible for the RR is the acceleration. In natural 
units, the magnitude of acceleration [m] and force [rri^] 
obtained by hitting a resting electron with a laser beam is 
usually not approaching unity (measured in units of the 
electron's mass). For example, for aq = 100 the maximal 
Lorentz invariant force is only \J —vP^Ua = 8 x 10~^ mo- 

However, in a head-on collision between the electron 
and the laser pulse, the "critical" force (i.e. acceleration 
of the order of the electron mass) can be achieved easily, 
as is presented in Fig. [31 This figure shows an acceler- 
ation of order unity, achieved by colliding a relativistic 



electron (70 = 1000) with a laser with intensity oq — 100. 
In fact, by computing the LF we see that "unity" accel- 
eration is achieved when 



a;ao7o ^ m^-, 



(48) 



corresponding to the appearance of the "critical" field 
strength E^ = m}(? jefi in the electron rest frame. This 
defines the solid (green) line in Fig. [2] Quantum effects 
may become relevant beyond this boundary [2ll |22| , and 
we have chosen examples within the domain where clas- 
sical dynamics certainly dominate. It is noteworthy that 
Eq. PSJ) also loosely demarcates the domain of validity 
of the LL equation itself, i.e., the lighter areas of Fig. [2] 
where the predicted radiation remains "small." For par- 
ticles other than the electron, condition Eq. (|15)) differs, 
as the constant oq defined in Eq. ([T5|) depends on the 
mass of the particle. 

Figure[5]exhibits the existence of an area of (ao, 70) pa- 
rameter space especially sensitive for experiments prob- 
ing radiation reaction. Between the solid (green) and 
dashed (blue) lines, the classical LAD and LL equa- 
tions are both valid, but their dynamical predictions 
in general can differ observably. The importance of 
this recognition is the reason we have chosen to provide 
in Sec. |V] numerical solutions for the point marked at 
(ao = 100,70 = 1000). 



B. Angular Distribution of Radiation 

The energy flux measured in the laboratory frame is 
2;iven by the Poynting vector [l[ 
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(49) 



where E', B are the electric and magnetic field respec- 
tively. Let y" = (t, y) be the spacctime point at which 
we evaluate the electromagnetic field, and denote by tq 
the proper time at which the worldline of the charge x" 
intersects the past light cone emanating from y". The 
retarded time tq is defined by the light-cone constraint 



f-a;o(To) = -R, 



(50) 



where we denote by R the spatial separation between the 
two events, namely 



i?= |j;-f(ro)|. 



(51) 



Since the Poynting vector in Eq. (^^ is quadratic in 
the fields, and radiation is defined as the energy trans- 
ported to infinity, only fields that are falling off as R~^ 
or slower radiate. Consequently, it can be shown [l| that 
the energy flux radiated to infinity by a single charge is 
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(52) 



where n is a unit vector in the direction of y — a:(r), 
w, a are the three-dimensional velocity and acceleration 
of the charge respectively, and subscript ret means that 
the quantity in the brackets is to be evaluated at the 
retarded time tq. 

The total amount of energy per unit solid angle emit- 
ted during a flnite acceleration period T is obtained by 
integrating the Poynting vector from the retarded lab 
time i = to t = r and is given by [l[ 
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(53) 



Since the solution in Eq. (1311) is given in terms of ^, it 
is useful to express the integrand using the variable ^ as 
wefl, as 
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In order to compute Eq. ([54|) we choose the coordinate 
system such that the electromagnetic wave propagates in 
the positive z direction while the particle initially moves 
in the negative z direction. Equation ([Ml) results in a 
cumbersome expression in the general case, which does 
not add much to the understanding of the physics of the 
problem. Instead, in the next section we will consider the 
angular distribution of the energy emitted in the speciflc 
cases presented earlier, of linearly and circularly polar- 
ized plane waves. 

Figure 2] shows the normalized angular distribution of 
radiation in the x — z plane in the laboratory frame for a 
LP plane wave hitting an electron initially at rest. The 
distribution was normalized to unity by dividing it by 
its maximal value in order to demonstrate the different 
shapes of radiation distribution for a laser beam with 
wavelength A = 942 nm and intensity oq = 0.1, 1, 10. For 
these intensities (and an electron at rest) the LF equation 
and the LL equation produce indistinguishable radiation 
distribution. The asymmetry in the dashed blue curve 
for ag = 10 is because the plot was generated for a laser 
wave duration of 26.8 fs. For these values of the wave 
duration and laser intensity, the velocity of the electron 
in the x direction does not finish an integer number of 
periods and results in an asymmetric motion in the x 
direction and consequently in the angular distribution of 
radiation. 



V. EXAMPLES 

We now discuss our main examples of an electron trav- 
eling in linearly and circularly polarized plane waves. 
Assume that in the lab frame, the electron is initially 
traveling towards the electromagnetic wave to a head-on 




500 



FIG. 4; The angular distribution of radiation for a linearly 
polarized wave. This is the normalized radiation distribution 
for an electron initially at rest, after interacting with a laser 
with ao — 0.1, ao = 1 and ao = 10 plotted in solid green, 
dotted red and dashed blue lines respectively. This plot is 
identical for the LF and for the LL equation. 



photon-electron collision at time i = 0. For the figures 
and the numerical values, we take an electromagnetic 
wave with intensity oq = 100, wavelength A = 942 nm, 
frequency to = 2 fs~ and a pulse that runs for a duration 
of 26.8 fs. The RR effects are demonstrated graphically 
for an electron initially at rest 70 = 1 and a highly rela- 
tivistic electron with initial 70 = 1000. 



A. Linearly polarized plane w^ave 

The definitions given in Eq. (|13p give the structure in- 
tegrals Eqs. dSni) and ((26)) 
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X" = I 0, -(2-3cosC-3Csin^ + cos''C),0,0 
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The projection of the four-velocity onto the wave vector 
is given by Eq. (fT9|) and yields 
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In Appendix 1X1 we provide the explicit expression for the 
four-velocity of the particle in this case, as is given by 
Eq.dSI]). 
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FIG. 5: The electron's Lorentz factor 7 and the longitudinal 
velocity Vz for a linearly polarized wave with ao = 100 and 
wavelength A = 942 nm hitting an electron initially at rest. 
The solid red line represents the solution of the LF equation, 
while the dashed blue line represents the LL equation which 
is indistinguishable in this case. 



Equation ([2T|) gives the relationship between ^ and r, 
which is in this case 

t(C) = -^ + iagro [1 + 2e ~ cos 2^] . (57) 

k ■ uq 8 

Notice that the limit tq — > indeed exists, where we 
obtain the linear relation of Eq. (|33l) . In fact, in this 
particular case of the LF equation with a linearly polar- 
ized plane wave, Eq. ([33)) is a known result [l^l, which 
means that the four-acceleration of the particle is always 
orthogonal to the wave vector k" in the absence of RR 
effects. 

Figure [5] shows 7 and the longitudinal velocity for an 
electron that is initially at rest. All figures are presented 
in the laboratory frame with respect to the laboratory 
time t. As Fig. [2] shows, this is still far from the RR 
dominated regime and the predictions of the LL equa- 
tion are practically indistinguishable from the ones given 
by the LF equation. The particle is boosted rapidly in 
the direction of the laser wave z and is rclativistic after 
approximately 100 as. 

If the particle is initially moving in the negative z di- 
rection, there is a head-on collision between the wave 
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FIG. 6: The electron's Lorentz factor 7, the transverse veloc- 
ity Vx and the longitudinal velocity v^ in the laboratory frame 
as a function of lab time for a linearly polarized wave with 
ao = 100 and an electron with initial 70 = 1000. The solid 
red line represents the solution of the LF equation, while the 
blue (dashed) line represents the LL equation. 



and the particle. An electron with initial 70 = 1000 
{Eq = 0.511 GcV) is on the critical line of the RR domi- 
nated regime. It is presented in Fig. |6l The pulse dura- 
tion is not long enough to flip the direction of the highly 
relativistic electron (a 180 fs pulse duration is needed). 
The radiation emitted by the particle slows the electron 
down rapidly in the direction of the wave, while it gains 
momentum in the direction of polarization. This is only 
true in the case of the LL model (the dashed blue line) 
due to the RR effects. In the absence of RR (the solid red 
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FIG. 7: The dependence between the laboratory time i and 
the variable ^ in a LP plane wave for LF (solid red line) and 
LL (dashed blue line) dynamics. The duration of ^ is doubled 
from the original duration of the wave in Lot. 
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FIG. 8: Radiation emission for a linearly polarized wave with 
ao = 100 colliding head-on with an electron with initial 70 = 
1000. The angle 6 is measured on the x—z plane, starting from 
the negative z axis. The radiation for the LF equation (solid 
red line) is 1 order of magnitude greater than the radiation 
in the case of the LL equation (the dashed blue line) . 



line) the pulse is incapable of slowing the particle down. 

Although the wave is active for about 8 periods, we 
see from Fig. [5] that the transverse velocity has 16 peri- 
ods, while the longitudinal velocity has 32 periods. The 
reason for this effect is twofold. First, recall that the so- 
lution Eq. (PT|) was given in terms of the variable ^, which 
has a nontrivial relation to the laboratory time in which 
we measure the velocities of the particle. This relation is 
demonstrated in Fig. [7] in which we see a doubling of the 
duration for ^. Second, a careful examination of the solu- 
tion for a linearly polarized plane wave (in Appendix \K\ 
shows that the period is doubled for the longitudinal ve- 
locity once more, resulting in 32 periods. The second 
effect is a feature of the linear polarization and is absent 
in the case of a circularly polarized wave. 

Since the electron is highly relativistic, the radiation 
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is emitted almost purely in the direction of its initial 
motion, namely, in the negative z direction. The dif- 
ferent radiation emission from the LF equation and the 
LL equation is shown in Fig. |5] where we see that the LF 
dynamics produce different radiation emission than the 
LL dynamics. 



B. Circularly polarized plane wave 

We use Eq. P^ to evaluate the structure integrals 
Eqs. dini) and ([26]) 



^ = -i 

X" = (O,sin(0 



Ccos(0,l-cos(0 



(58) 
Csin(C),0). 



The projection of the four-velocity onto the wave vector 
is given by Eq. (fTQ]) and yields 



k ■ u = 



k ■ uq 



1 + Toa5(fc • uo)^ 



(59) 



The four- velocity of the particle is given by Eq. ([5T|) and 
the analytic expression for it in this case is given in Ap- 
pendix |X1 In order to have the four- velocity in terms 
of the proper time r, wc need to invert Eq. (j2ip . In the 
case of a circularly polarized plane wave, we can find the 
inverse explicitly and obtain 



e(r) = 



^1 + 2a§(fc ■ Mo)Vor - 
al{k ■ uo)to 



(60) 



Since a circularly polarized wave with intensity Oq = 
100 produces similar results to a linearly polarized plane 
wave of the same intensity (Fig. [S]), the graphs are not 
presented here. However, because for the circularly po- 
larized case in Eq. ([T?]) there is no choice of phase ^o of 
the wave for which the transverse motion of the resting 
electron vanishes, the phase at the instant of collision has 
physical significance resulting in a drift in the direction 
of initial polarization (x in this case) . This drift is a non- 
physical consequence of the nonphysical electromagnetic 
field. A different choice of initial phase for the linearly 
polarized wave produces a drift in the direction of polar- 
ization as well. This mathematical artifact is a result of 
a plane wave whose duration is infinitely long, as one can 
see from Eq. (|3T|) that the solution depends explicitly on 
the four-potential at the initial time Aq. If the initial 
four-potential Aq is nonzero, the second term in Eq. ([31]) 
causes a drift in the opposite direction of Aq . 

In practice, any experiment involving laser beams will 
require a certain finite duration on which the wave is 
being turned on. This can be modeled into the solution 
by defining a turn-on function 
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which vanishes for negative ^, and increases continuously 
to 1 over a finite period determined by the parameter a. 
The new shape function for the four- vector potential is 
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for the circularly polarized wave [cf. with Eq P^ ]. The 
introduction of the turn-on function H{^) in the shape 
function guarantees that the electromagnetic field will be 
turned on smoothly. The integration appears since /(^) 
defines the four-potential, but we would like a smooth 
turn-on of the electromagnetic field which is defined in 
terms of the derivative of /(^) [see Eqs. ^ and ([TT]) ]. 

In the following graphs we choose the parameter a to 
be such that the electromagnetic field is turned on on 
a period of about a quarter of a wavelength ^ A/4 == 
210 nm. Since the turn-on function makes the mathemat- 
ical expressions in the solution of Sec. IIIII cumbersome, 
they were omitted. 

Figure IHl presents the head-on collision between the 
wave and a relativistic particle (70 = 1000). The circular 
polarization is much more efiicient in slowing the parti- 
cle down, as can be seen by comparing the velocity in 
the z direction in the two cases (compare Figs. Inland [5]). 
About 26 fs suffice for a circularly polarized wave to stop 
the highly relativistic electron, after which the electron 
reverses its direction of motion to coincide with the direc- 
tion of propagation of the wave. Figure [TU] presents the 
radiation emission, where we see that the LF equation 
results in the familiar dipole radiation (the angle of sep- 
aration between the two branches of the dipole is of order 
1/7 ^ 10""^). Similarly to the LP case, the LL equation 
produces radiation which is 1 order of magnitude smaller 
than the LF equation. 



VI. CONCLUSIONS 

The dynamics of a radiating charged particle has 
been the subject of investigation for more than a cen- 
tury. Lorcntz himself recognized the need to expand 
the MaxwcU-Lorentz framework to include an additional 
force due to radiation reaction. In the intervening years, 
many others sought physical consistency by including a 
radiation-reaction force as a modification of the Lorentz 
force equation. These efforts have produced numer- 
ous radiation-reaction models but unlike the Maxwell- 
Lorentz equations none of these has been derived from 
an action principle. 

So far. the theoretical efforts have remained experi- 
mentally untested. Our objective here has been to con- 
nect theory with possible experiment. We found that 
when colliding relativistic electrons (7 > 1000) with an 
extremely intense laser beam (oq > 100, power > 2 x 10^^ 
W/cm^), the laser is capable of stopping the relativistic 
electron once radiation-reaction force is included. 

For the linearly polarized wave, the chosen pulse length 
was insufficient to reverse the initial momentum of the 
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FIG. 9: The electron's Lorentz factor 7, the transverse veloc- 
ity Vx and the longitudinal velocity v^ for a circularly polar- 
ized wave with ao = 100 hitting an electron with 70 — 1000. 
The solid red line is the solution of the LF equation, while 
the dashed blue line is the LL equation. 



electron. However, a longer laser pulse could reverse the 
net momentum of the electron, accelerating it back in the 
direction from which it came. Thus being a dissipative 
effect, radiation reaction makes the gross consequences 
sensitive to the duration of the interaction Q. A more 
thorough exploration of the outcomes and observables for 
various pulse lengths using the solutions here presented 
would be an important preliminary to any experiment. 

We have presented an analytic solution of the Landau- 
Lifshitz Eq. ([5]) , and an analysis of the radiation distri- 
bution emitted by the charged particle. The radiation- 
reaction- dominated regime criterion given by Eq. (j4ip is 
visualized in Fig. [21 showing that the radiation-reaction 
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FIG. 10: Radiation emission for a circularly polarized wave 
with ao = 100 colliding head-on with an electron with initial 
7o = 1000. The angle 9 is measured on the x — z plane, 
starting from the negative z axis. The radiation for the LF 
equation (solid red line) is 1 order of magnitude greater than 
the radiation in the case of the LL equation (the dashed blue 
line). 



force can be probed using current pulsed laser systems. 
Figure [H exhibits the domain in which the classical dy- 
namics here solved dominates, but the Landau-Lifshitz 
equation results in different radiation patterns than those 
obtained from the Lorentz force Eq. ^. This suggests 
that a study of the radiation patterns emitted in such 
experiments will provide a first experimental opportu- 
nity, well within the reach of current laser systems, in 
which radiation-reaction effects can be observed and un- 
derstood in the new physics domain. 

Examining the radiation emission in both Figs. [8] and 
[TOl one sees that the LL equation predicts less total in- 
tegrated radiation emission than the LF equation. On 
first sight this appears to conflict with rapid, large am- 
plitude oscillations of the electron. However, we recall 
that the radiation rate is proportional to the accelera- 
tion, TZ ^ ii^, and Fig. [3] for \J —v? or Figs. [HI and IH] for 
Lorentz factors exhibit how the LL equation permits the 
electron to decelerate while the LF does not. One may 
thus expect the energy converted to radiation by the ac- 
celerated motion of the electron to decrease with time 
for the LL equation, but to remain constant for the LF 
equation. As a result, the total integrated emitted power 
predicted by the Lorentz force can indeed be larger. For 
this reason, one cannot evaluate the radiation produced 
in high-intensity laser-matter interactions according to 
the Lorentz force dynamics alone. 

The examples here thus exhibit the more general point 
that accurate predictions of radiation emission require in- 
corporating radiation-reaction effects whenever high ac- 
celerations are expected. It should be emphasized that 
our radiation results are computed according to Dirac's 
evaluation of the radiated momentum, which leads in 
turn to the LAD dynamics. On this point, an im- 
provement on the present work would be derivation of 
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a Landau-Lifshitz-like equation for particle motion self- 
consistent with the particle's radiation field. 

The authors hope that this work will assist in explor- 
ing the nature of radiation-reaction effects on firm experi- 
mental ground, and support the search for a fundamental 
theory of electromagnetism applicable in the high accel- 
eration regime. In that way the investigations presented 
here also reveal the possibility of reaching for fundamen- 
tal physics horizons. When the applied force approaches 
unity in natural units [m?], the current physics frame- 
work has not been tested and we are at the limit of our 
understanding. 
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Appendix A: Supplemental formulas for the 
examples 
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a) Linearly polarized plane wave 



Equation pip for the linearly polarized plane wave in 
Eq. ^ gives 
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b) Circularly polarized plane wave 

The solution in this case is 
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ao(l-cosO + ao (2-2cos^-^sinC+ y) +ao (1 -cos^-^sin^ 



f 



^2M 

2 



(A2) 
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Appendix B: The four-acceleration 

We differentiate Eq. (|3T|) with respect to t [remember that in Eq. ((3T|) the four-velocity is given in terms of ^] and 
obtain the four-acceleration of the particle according to the LL equation: 



k ■ u 



k ■ Mo 

+To{k-uf 



{(fc.Mo)aoi'" 
aoi"" 



aoi'■wo + ao(^-^)•^']fc"} 



2 ^ ■ " A'^u^ 



k ■ u 



A'"(A"-AE^)-a;5V^" 



+ 



aoi" • uo - alA'^ - alA' ■ {A' - if,) - al{A - Aq) ■ A" + a^^/jA' ■ uq 



3 ^ ■ " A'^iA - Ao) ■ uo + 4A' ■ X + 4iA - io) • i'V- - a^^i'^ ^^ "^"^^ 



fc • Uq 

+T^{k ■ uf\ \al{k- M)i'2(i'" - i(f ) - a^(fc • u)i'2x" 



A: • Uq 



k" 

k ■ Uq 



al{A'-A',)-A" 



k ■ u 



k ■ Uq 



A'\A' - 4) • UQ + 4A" ■ X + 4{A' - 4) • iV + 4^A'^ - a4^ji^,2^ 



k ■ Uq 



k- 



k ■ u 



k ■ u 



-aQ- A'\A - Aq) ■ {A' - A'q) + al- A'\ ■ uq - a^x ■ A'^ + a^- A'^{A ~ Aq) ■ x 

k ■ Uq k ■ Uq k ■ Uq 



k" 



3 Ai2 



+T^ik-uyA 



at^-^^^+al{A'-A'Q)-X + a' 






«^ Y ^fc". 



(Bl) 



where we collected terms in powers of tq and separated 
them according to the direction of propagation and de- 
pendence on the intensity of the wave Uq as before. 



Appendix C: Comparison with Previous Works 
1. Analytic Comparison 

As an initial value problem, the Landau-Lifshitz 
Eq. (O satisfies the conditions of the uniqueness and ex- 
istence theorem. Therefore, one expects the evolution of 
the particle to be uniquely determined by the equation 
together with the initial conditions. This raises the ques- 
tion as to whether the solution given here by Eq. pip and 
the solution given in |16j are equivalent or even consis- 
tent. 

The differences in notation between our paper and [l6| 
are presented in Table IC II and can be used to identify 
our expressions with the equivalent ones given in jl6| . 



TABLE I: Notations in our paper compared with [Tg| 



This paper 


Reference [16] 


k" 


n" 


^ = fc • u 


(j) — n ■ u 


k ■ uo 


Po 


k-u 
k-UQ 


h(4>) 



In order to see equivalence between the two solutions 
wc define the four-vector 

I"(e) = aoAA"(e) + (fc.Mo)ToaoA^'"(C) (CI) 
-(fc.Mo)Toagx"(C) 
and the scalar quantity 



Ki) = ^ = 1 - (fc . UQ)TQali:{0- 
They allow us to express the solution in Eq. pip as 



(C2) 



k ■ u 
k ■ Uq 

1 



2k ■ Uq 



l]fc" 



(C3) 



k ■ Uq 



[k"I^ ~ fc'^Z"] Uq, 
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2k ■ Uq 



-I^k" 



where I" and h are functions of the phase £,. Using Table 
IC 1\ this can be compared with the solution given in Eq. 
(ll)of [l|: 

."^i{< + ^[/,2-l]n<^ (C4) 



1 

Po 



ra^ 



i,^i^+i^^'^ 



eff 



U0,/3 



-^[«;if + gi|]„-}, 



where 

1(0) 



dkp 



■' 6 m-' 



(C5) 
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and 



/•a/3 



(C6) 



Here a" are two four-vectors that can be associated with 
the polarization of the wave, and n" is the wave four- 
vector. 

To show the equivalence of Eq. (|C3|) and Eq. ()C4|) . we 
rewrite four- vector ()Cip using Eq. ([26l) as 



Z"(0 



aoAi"(0-(fc-wo)Toai]x"(0 (C7) 

+(fc-7/o)ToaoAi'"(0 

ao / A'"{y)[l~{k-u^)Toalij{y)\dy 
Jo 







ao 



A"'{y)h{y) + {k-uo)T^A"{y) 



dy. 



This reveals that I" in our notation is nothing other than 
Xia" +X2a2 in the notations of [la |. 

Therefore although the notations in this paper and 
paper [lal differ, Table IC II together with the result of 
Eq. (|C7p show that the two solutions in Eq. (|C3|) and 
Eq. (|C4p are identical and both solve the Landau-Lifshitz 
Eq.®. 



2. Numerical Comparison 



In 27[ the influence of radiation reaction according to 
the LL equation was studied numerically for a head-on 
collision of an electron with a strong laser pulse. The 
electron was initially traveling in the negative y direc- 
tion with initial energy of 40 MeV and the laser wave 
propagating in the positive y direction with frequency 
u) = 1.55 eV and an intensity /g = 5 x 10^^ W/cm^. The 
laser wave was linearly polarized in the z direction and 
subject to a sine-squared envelope with duration of 27 fs, 
resulting in a motion in the y — z plane given by Fig. 1 

of [13. 

In our notation, the above parameters translate to w = 
2.35 fs"\ ao = 152.925 and 70 = 78.28. The envelope can 
be modeled into the analytic solution (|5T|) by defining [see 
Eq. dH)] 



envelope shape (the factor 20 was found by comparing it 
to the image in [231, ^^ the exact value was not stated). 
Figure [11] shows the corresponding analytical result ob- 
tained from Eq. pip , and shows the same qualitative 
features as pT\ . 

The trajectories in [23| and Fig. [TT] cannot match per- 
fectly, as [27| in principle considers a focused laser beam, 
which would have a nonplanar wave front that is not 
implemented in our analytic solution. Notably, we do 
not find any transverse drift comparable to Fig. 1 of 
[27| . Nevertheless, the solution provided here is versa- 
tile enough to be used to provide the same features that 
one would expect from a focused beam. This capability 
is a result of the containment apparent in Fig.[TT] of the 
transverse oscillation to < 0.2A from the beam axis, a 
distance over which the envelope of a focused laser pulse 
does not appreciably vary. 



z[A] 




y[A] 



-3.5 -3.0 -2.,5/ 



iO;:;;:-rl;.5: 



V.O '--^0.5 



0.15 : 

0.10 ': 
0.05 ': 

-0.05 ': 
-0.10 ': 

-0.15 '-- 



yW 



/(O 



cos{y)ain^{y/20)dy, 



(C8) 



where the cosine in the integral gives the usual linear po- 
larization, while the sine squared term gives the required 



FIG. 11: The analytical results for ao = 152.925 and 70 — 
78.28 for LF (top panel) and LL (bottom panel) solutions. 
Compare to the numerical result in Fig. 2 of [27| which demon- 
strates similar qualitative features. 
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